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o : 

, We consider a system consisting of an atom in the approximation of a harmonic oscillator of 

frequency Oj, coupled to the scalar potential inside a spherical reflecting cavity of radius R. We 
use dressed states introduced in a previous publication [Andion, Malbouisson, and Mattos Neto, J. 
, Phys. A 34, 3735 (2001)], which allow a non-perturbative unified description of the atom radiation 

' process, in both cases, of a finite or an arbitrarily large cavity. We perform a study of the energy 

^SJ ' distribution in a small cavity, with the initial condition that the atom is in the first excited state 

and we conclude for the quasi-stability of the excited atom. For instance, for a frequency ub of the 
order lo ~ 4.00 x 10^*/s (in the visible red), starting from the initial condition that the atom is in 
the first excited level, we find that for a cavity with diameter 27? ~ 1.0 x 10~^m, the probability 
that the atom be at any time still in the first excited level, will be of the order of 97%. For a typical 
microwave frequency ~ 2, 00 x 10^''/s we find stability in the first excited state also of the order 
of 97% for a cavity radius ii ~ 1.4 x 10~^m. 
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I. INTRODUCTION 



C/2 

' There are situations in the domain of Atomic Physics, Cavity Electrodynamics and Quantum Optics, where per- 
O • turbation methods are of little usefulness, for instance, in resonant effects associated to the coupling of atoms with 
strong radiofrequency fields The theoretical understanding of these effects on perturbative grounds requires the 
calculation of very high-order terms in perturbation series, what makes the standard Feynman diagrams technique 
practically unreliable in those cases Ipj. The trials of treating non-perturbativelly such kind of systems, have lead to 



. the idea of dressed atom, introduced in refs and Since then this concept has been used to investigate several 
■ situations involving the interaction of atoms and electromagnetic fields ( Q], ||], [^). A way to circumvect the mathe- 
. matical difficulties due to non-linear character of the problem, is to assume that under certain conditions the coupled 
^-H . atom-electromagnetic field system may be approximated by the system composed of an harmonic oscillator coupled 
^""^ ' linearly to the field trough some effective coupling constant g. This is the case in the context of the general QED 
i linear response theory, where the electric dipole interaction gives the leading contribution to the radiation process ( 
^ 0, §, @) and also in several branches of Quantum Optics ( |ll|, [Q, Q). 

^ . We consider a system composed of an atom (approximated by a harmonic oscillator) coupled linearly to the scalar 
. ^ ' potential, the whole system being confined inside a reflecting sphere of radius R. We give a non-perturbative treatment 
to the field-atom system introducing some dressed coordinates that allow to divide the coupled system into two parts, 
the dressed atom and the dressed field, what makes unnecessary to work directly with the concepts of bare atom, bare 
^ field and interaction between them. For instance, to describe the radiation process, having as initial condition that 
• • only the mechanical oscillator (the atom), be excited, the usual procedure is to consider the interaction term in 
. 5^ ' the Hamiltonian written in terms of go and the field modes qi as a perturbation, which induces transitions among the 
, eigenstates of the free Hamiltonian. In this way it is possible to treat approximately the problem having as initial 
\^ ' condition that only the bare mechanical oscillator (the atom) be excited. But as is well known this initial condition is 
5^ , physically not consistent in reason of the divergence of the bare oscillator frequency, due to the interaction with the 
field. The traditional way to circumvect this difficulty is by the renormalization procedure, introducing perturbativelly 
order by order corrections to the oscillator frequency. In this paper we adopt an alternative procedure introduced in 
. We do not make explicit use of the concepts of interacting bare oscillator and field, described by the coordinates 
go and {qi\. We introduce dressed coordinates q^ and {g^} for, respectively the dressed atom and the dressed field 
modes. In terms of these new coordinates a non-perturbative approach of the radiation process and of the distribution 
of energy inside the cavity is possible. 
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We start considering an atom approximated by a harmonic oscillator qq (t) of frequency ujq (we will introduce below 
a renormalized frequency uj which is physically meaningfuU) coupled linearly to the scalar potential 0, the whole 
system being confined in a sphere of radius R centered at the oscillator position. The equations of motion are, 



qo(t) + ujlqoit) = 2TTy^ / d'^r(f>{r, t)S{r) 



(1.1) 



(1.2) 



Using a basis of spherically symmetric Bessel functions defined in the domain < |r| < i?, the equations above 
can be written as a set of equations coupling the atom to the harmonic field modes, which can be derived from the 
Hamiltonian 
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(1.3) 



II. THE EIGENFREQUENCIES SPECTRUM 

We consider for a moment as in |lj], the problem of a harmonic oscillator go coupled to N other oscillators. In the 
limit N oo we recover our original situation of the coupling oscillator-field after redefinition of divergent quantities, 
in a manner analogous as renormalization is done in field theories. In the above equations, g is a coupling constant 
(with dimension of frequency), rj = y/2gSuj and Alu = ttc/R is the interval between two neighbouri ng fi eld frequencies, 
LUi^i — tjji = Auj = -Kc/ R and qi stands for the harmonic modes of the field. The Hamiltonian (1.3) can be turned 



to principal axis by means of a point transformation, = t^Qr , P/i = ^^-fr, performed by an orthonormal matrix 
T = (ip, ^ = (0,fc), k — 1,2, ...iV, r = 0,...A^. The subscript and k refer respectively to the atom and the 
harmonic modes of the field and r refers to the normal modes. The transformed Hamiltonian in principal axis reads, 
H = \ Er=oiP? + ^rQl), where the fir's are the normal frequencies corresponding to the possible collective oscillation 
modes of the coupled system. The eigenfrequencies fir satisfy the equation , 



-'.-Nv'-n^ = v'f:^^- (2-1) 



The + 1 solutions fir of Eq.(2.1), correspond to th e + 1 normal collective oscillation modes. 

It can be shown flj] that if cjq >Nrf', Eq.(2.1) yields only positive solutions for fl^ (all collective modes are 
harmonic), while if < Nif ^ Eq.(2J^) has a negative solution fi^_. This means that there is a damped oscillation 
mode that does not allows stationary configurations. We will not care about this last situation. Nevertheless it should 
be remarked that in a different context, it is precisely this negative squared frequency solution (runaway solution) 
that is related to the existence of a bound state in the Lee-Friedrechs model. This solution is considered in ref. p^ ] 
in the framework of a model to describe qualitatively the existence of bound states in particle physics. Thus we take 
uJq > Nr]^ and define the renormalized oscillator frequency lj, Co — ■\/wq — Nrf . In the limit N ^ oo the meaning 
of the frequency renormalization becomes clear. It is exactly the analogous of a mass renormalization in field theory, 
the infinite too being chosen in such a way as to make the renormalized frequency Cu finite and equal to the observed 



oscillator frequency. In terms of the renormalized frequency Eq.(2.1) can be writen, after some manipulations, in the 
form [l^ , 

,Rn, fl c , 

cot — = — + — 1 . 2.2 

C TTg Rii TTgC 

The solutions of Eq. ( ^.^ ) with respect to fl give the spectrum of eigenfrequencies fir corresponding to the collective 
normal modes. The transformation matrix elements turning the oscillator-field system to principal axis is obtained 
taking the limit N —> oo, after some rather long but straightforward manipulations in |p^. They read. 
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We define below some coordinates ^q, q[ associated to the dressed atom and the dressed field. These coordinates 
will reveal themselves to be appropriate to give an appealing non-perturbative description of the atom-field system. 



III. DRESSED STATES 



We start from the eigenstates of our system, represented by the normalized eigenfunctions, 



^l«....(0,^)-^ 



Toe" 



(3.1) 



where stands for the n^-th Hermite polynomial and Tq is the normalized vacuum eigenfunction. Let us introduce 
dressed coordinates Qq and {q[} for, respectively the dressed atom and the dressed field, defined by jlj]. 



(3.2) 



valid for arbitrary R and where = {w, toi}. In terms of the bare coordinates the dressed coordinates are expressed 
as. 



(3.3) 



Let us define for a fixed instant the complete orthonormal set of functions 



■^pKoK^..iq') = n 




ro, 



(3.4) 



where q'^ — (jg, q^, oj^ = {uj, uji\. Note that the ground state Fq in the above equation is the same as in Eq.( [3.l[ ) 
The invariance of the ground state is due to our definition of dressed coordinates given by Eq.(3.2). Each function 
...il') de scribes a state in which the dressed oscillator q'^ is in its — t h ex cited state. Using Eq.(3.2) the 



functions (3.4) can be expressed in terms of the normal coordinates Qr- But since (3.1) is a complete set of orthonormal 
functions, the functions (3.4) may be written as linear combinations of the eigenfunctions of the coupled system (we 
take t = for the moment), 

V'«o«....('?')= E T^"o""r.-(o)'^"o«,«....(Q,o), (3.5) 

no"!--- 

where the coefficients are given by, 

TZZ\::{o) = J dQ^„„«,...0„„„,„,... , (3.6) 

the integral extending over the whole Q-space. 

We consider the particular configuration ip in which only one dressed oscillator q'^ is in its iV-th excited state, all 
other being in the ground state. 



^o...NMo.M) = (2-^iV!)-^i/^(^/^^;)ro 



The coefficients ( p.6[ ) have been calculated in ref. ||lj]. We get, 

, Nl 
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(3.7) 



(3.8) 



where the subscripts /i = 0, i refer respectively to the dressed atom and the harmonic modes of the field and the 
quantum numbers satisfy the constraint no + ni + ... — N. In the following we focus our attention on the behaviour of 
the system with the initial condition that only one dressed oscillator q'^ (the dressed atom or one of the modes of the 



dressed field) be in the iV-th excited state. We will study in detail the particular case N = 1, which will be enough to 
have a clear understanding of our approach. Let us call t he coiifiguration in which only the dressed oscillator q'^ is 



in the first excited level. We have from Eqs.(p.7| 
of the first-level excited dressed oscillator q'^, 



(3.5) (3.8) and (3^) the following expression for the time evolution 



(3.9) 



where the coefficients f^'^{t) are given by 



(3.10) 



From Eq.( |3.9[ ) we see that the initially excited dressed oscillator naturally distribut es its energy among itself and all 
other dressed oscillators as time goes on, with probability amplitudes given by Eq.(3.10). If the dressed oscillator q^ 
(the atom) is in its first excited state at t = 0, its decay rate may evaluated from the time evolution equation, 



r;(t)-E/"''Wri(o) 



(3.11) 



In Eq.(3.11) the coefficients f^^{t) have a simple interpretation: f^^{t) and f^^t) are respectively the probability 
amplitudes that at time t the dressed atom still be excited or have radiated a photon of frequency uji. We see that this 
formalism allows a quite natural description of the radiation process as a simple exact time evolution of the system. 
We consider in the following the time evolution of the excited atom, in the cases of a very large and a very small 
cavity. 

A very large cavity: 

In the case of a very large cavity our method generalizes what can be obtained from perturbation theory. The 
probability that the atom be still excited at time t can be obtained in continuous language from the amplitude given 
by Eq. (Pi, 
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(3.12) 



For large t {t >> ^), but for in principle arbitrary coupling g, we obtain for the probability of finding the atom still 
excited at time t, the result fl^ . 
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(3.13) 



p-. In the above expression the approximation t » ^ plays a role only in the two last terms. 



where uo — 

due to difhculties to evaluate exactly the integral in Eq. ( ^.12 ) along the imaginary axis using Cauchy's theorem 



The first term comes from the residue at VL — lj -\- and would be the same if we have done an exact calculation. 



If we consider in eq. (3.13) g « cD, which corresponds in electromagnetic theory to the fact that the fine structure 
constant a is small compared to unity (for explicit calculations we take below g/u) = a), we obtain the well known 
perturbative exponential decay law. 



IV. THE RADIATION PROCESS IN A SMALL CAVITY 



Let us now consider the atom placed at the center of a very small cavity, i.e. that satisfies the condition that its 
radius be much smaller than the coh eren ce lenght, R << c/g. To obtain the eigenfrequencies spectrum, we remark 
that from a numerical analysis of Eq. (|2.2| ) it can be seen that in the case of a small cavity radius i?, its solutions are 
near the frequency values corresponding to the asymptots of the curve cot(^), which correspond to the field modes 
Ui = iiTc/R. The smallest solution departs more from the first asymptot than the other larger solutions depart from 
their respective nearest asymptot. As we take larger and larger solutions, they are nearer and nearer to the values 
corresponding to the asymptots. For instance, for a cavity radius R of the order of 10~^m and oj ^ 10^°/s, only the 
lowest eigenfrequency fio is signicantly different from the field frequency corresponding to the first asymptot, all the 



other eigenfrequcncies ^Iki k — 1,2, ... being very close to the field modes knc/R. For higher values of ui (and lower 
values of R) the diffe renc es between the eigenfrequcncies and the field modes frequencies are still smaller. 

Thus to solve Eq.(2^) for the larger eigenfrequcncies we expand the function cot(^) around the values corre- 
sponding to the asymptots. We write, 



TTC 



k^l,2,.. 



with < ek < 1, satisfying the equation, 



cot(7rek) = — (k + Ek) 
gR 



1 



(k + ek) 



TTgC 



But since for a small cavity every ek is much smaller than 1, Eq.(4.2) may be linearized in ek, giving, 

TTgcRk 



(4.1) 



(4.2) 



(4.3) 



Eqs.(ll) and (13) give approximate solutions to the eigenfrequcncies flk, k = 1,2.... 

To solve Eq.(2.2) with respect to the lowest eigenfrequency f^o, let us assume that it satisfies the condition Q.aR/c « 
1 (we will see below that this condition is compatible with the condition of a small cavity as defined above) . Inserting 
the condition Q.qR/ c « \ in Eq.( ^.2| )and keeping up to quadratic terms in the solution for the lowest eigenfrequency 
r^o can be writen, 



(4.4) 



Consistency between Eq.(4.4) and the condition QqR/c << 1 gives a condition on the cavity radius, 
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(4.5) 



Let us define the coupling constant g to be such that g = Qa, where a is the fine structure constant, a — 1/137. Then 
the factor multiplying c/g in the above equation is ~ 0.07 and the condition R c/g is replaced by a more restrictive 
assumption R ^ 0.07{c/g). For a typical infrared frequency, for instance cZ; ~ 2,0 x 10^^/s, our calculations are valid 
for a radius R ^ 10^'^m. 

From Eq.(2.3) and using the above expressions for the eigenfrequcncies in a small cavity, we obtain the matrix 
elements. 
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(4.6) 



To obtain the above equations we have neglected the corrective term ek , from the expressions for the eigenfrequcncies 
rifc. Nevertheless, corrections in should be included in the expressions for the matrix elements i^, in order to avoid 

spurious singularities due to our approximation. 

Let us consider the situation where the dressed atom is initially in its first excited level. Then from Eq.( |3.1C| ) we 
obtain the probability that it will still be excited after a ellapsed time t. 



= {tlf + 2Y,{Af{A? cosiiik - n„)t + J2 (to fit'o)' cos(r!fe - not 
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Using Eqs.(O) in Eq.(4.7), we obtain 



\f°{t)\^ = l-nS + Ai- 



oo I 4^1 



fc=l 
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(4.7) 



(4.8) 



where we have introduced the adimensional paramet er S = R g/c <^ 1, corresponding to a small cavity and we 
remember that the eigenfrequcncies are given by Eqs.( [4.ll ) and (4.3). As time goes on, the probability that the atom 
be excited attains periodically a minimum value which has a lower bound given by. 



Min(r(t)n = l-^^+l|L!^^ (4.9) 
For a frequency ui of the orde r (D ~ 4.00 x 10^'*/s (in the red visible), which corresponds to S ^ 0.005 and 2R ~ 



1.0 X 10^ TO, we see from Eq.(O) that the probabihty that the atom be at any time excited will never fall below 
a value ~ 0.97, or a decay probability that is never higher that a value ~ 0.03. In other words, atoms having such 
emission frequency, placed in a such a small cavity in the first excited level, will be stable in the excited state to the 
order of 97%. It is interesting to compare this result with experimental observations in |]l6|, where stability is 
found for atoms emiting in the visible range placed between two parallel mirrors a distance L = 1.1 x 10~^to apart 
from one another. For lower frequencies the size of the cavity ensuring quasi-stability of the same order as above, for 
the excited atom may be considerably larger. For instance, for w in a typical microwave value, (D '--^ 2, 00 x 10^°/s and 
taking also S ~ 0.005, the probability that the atom remain in the first excited level at any time will be larger than 
a value of the order of 97%, for a cavity radius R ~ 1.0 x 10~^to. The probability that the atom remain excited as 
time goes on oscillates with time between a maximum and a minimum values and never departs significantly from the 
situation of stability of the atom in the excited state. Indeed for an emission frequency Co ~ 4.00 X lO^-^/s (in the red 



visible) considered above and i? ^ 1.0 x 10 to, the period of oscillation between the minimum and maximum values 

j_ 

12 

the period is T ^ ^ x 10"^°s. 



of the probability that the atom be excited, is T '--^ x 10 ^^s, while for w ~ 2, 00 x 10^°/s, and i? 1.4 x 10 ^to. 



V. CONCLUDING REMARKS 



We have used in this paper a formalism that allows an unified approach to the radiation process by an atom, in 
rather different situations, as the atom confined in a very small cavity or in free space. The behaviour of atoms 
confined in small cavities is completelly different from the behaviour of an atom in free space or in a large cavity. In 
the first case the emission process is very sensitive to the presence of boundaries, a fact that has been pointed out 
since a long time ago in the literature ( jl^ , jl^ , pof ) . Our dressed states approach gives an unified description for 
the dressing of the atom by the field modes and the emission process in a cavity of arbitrary size, which includes 
microcavities and very large cavities (free space emission). We recover here with our formalism the experimental 
observation that excited states of atoms in sufhciently small cavities are stable. We are able to give formulas for the 
probability of an atom to remain excited for an infinitely long time, provided it is placed in a cavity of appropriate 
size. For an emission frequency in the visible red, the size of such cavity is in good agreement with experimental 
observations ( |16|, Q). Also, our approach gives results in good agreement with previous theoretical results for the 
emission in free space, generalizing the well known exponential decay law. Moreover the detailed behaviours which we 
obtain with our formalism are very different in the two situations: The atom in a very large cavity has a probability 
decay rate weekly oscillating and monotonically varying with time (the probability that the atom be excited decreases 
almost exponentially with increasing time). In the case of an excited atom placed in the center of a very small cavity, 
the probability that it remains excited as time goes on oscillates very rapidly with time and never departs significantly 
from the situation of stability of the atom in the excited state. 
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